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Abstract 

The aim of this work is to determine the thermal conductivity of mesoporous silicon (PoSi) by fitting the 
experimental results with simulated ones. The electrothermal response (resistance versus applied current) of 
differently designed test lines integrated onto PoSi/silicon substrates and the bulk were compared to the 
simulations. The PoSi thermal conductivity was the single parameter used to fit the experimental results. The 
obtained thermal conductivity values were compared with those determined from Raman scattering 
measurements, and a good agreement between both methods was found. This methodology can be used to easily 
determine the thermal conductivity value for various porous silicon morphologies. 
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Background 

Application of porous silicon for electrical isolation in 
radiofrequency devices is widely studied [1,2]. However, 
characterization of its thermal properties is very import- 
ant, too, in order to anticipate any thermal risk due to 
self-heating effects. Thermal conductivity of porous sili- 
con depends mainly on its own morphological para- 
meters at nanoscale level, such as pore morphology and 
average size of Si nanocrystallites [3] . 

There are several existing methods to measure thermal 
conductivity which include scanning thermal micros- 
copy, that maps the local temperature and thermal con- 
ductivity of an interface [4], the 3o) method, which 
determines the thermal conductivity by applying an AC 
signal to metal strip [5], or the modeling approach [6], 
which is used in modeling the effect of pore size, pore 
arrangement, and porosity. The principle of comparing 
electrothermal simulations and experiments in order to 
determine thermal parameters of a material has already 
been applied successfully [7,8]. 



In this study, we propose to apply this method in 
order to evaluate the thermal conductivity of a porous 
silicon nanomaterial (mesoporous in our case) by elec- 
trical characterization of metal strips with subsequent 
electrothermal finite element simulations. These metal 
strips are integrated with the dielectric layers on top of 
bulk silicon and mesoporous silicon PoSi/silicon sub- 
strates. The electrical tests have been performed with a 
probe station on various test lines with different widths 
of 10, 30, 50 \im and a constant length of 1,100 (im. The 
electrical resistance of the test line is measured versus 
the applied current. This resistance is not constant with 
the current and varies due to the resistance temperature 
dependence. Firstly, the model parameters have been 
calibrated by comparing experimental and simulation 
results for test lines integrated onto bulk silicon sub- 
strate. Secondly, thermal conductivity of the mesoporous 
silicon has been determined by fitting the experimental 
and simulation results for test lines integrated onto 
PoSi/Si substrate. To validate this approach, a micro- 
Raman thermal conductivity measurement is performed. 
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Methods 

Silicon anodization 

Six inch porous silicon samples were fabricated by 
anodization of highly w-doped (100) silicon wafers 
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(p = 10 mH-cm). The anodization was run in a HF 
(50%) -acetic acid- water solution (4.63:2.14:1.43). Acetic 
acid has been used as a surfactant which improves the 
electrolyte penetration into the pores [9]. Moreover, acetic 
acid is known to improve the roughness and the uniform- 
ity of porous layers as well as to increase the etch rate [10]. 

Anodization was performed in a double tank electro- 
chemical cell. A mesoporous layer 100-|im thick was 
obtained by 60 min of anodization at 14.5 mA/cm 2 . 
Mesoporous thicknesses were measured by scanning 
electron microscopy (SEM) after splitting the samples. 
An average 20% porosity was evaluated by weight mea- 
surements. Finally, porous silicon samples have been 
annealed at 300°C under N 2 for 1 h in order to stabilize 
the structure. The typical morphology observed in our 
case is presented in Figure 1. 

Test line patterns 

Test lines with various designs were integrated on the 
porous silicon samples of highly w-doped silicon (p = 10 
mO-cm). A 500-nm plasma enhanced chemical vapor 
deposition (PECVD) oxide has been deposited as a cap 
layer on the porous silicon. Test lines were made with 
l-|im aluminum deposited at 350°C by physical vapor 
deposition (PVD). A 500-nm PECVD oxide was then 
deposited to passivate the aluminum. Finally, the elec- 
trical contact to the first metal level has been made with 
l-|im PVD aluminum. The patterns were defined by 
standard photolithography and dry etching with chloride- 
based plasma by reactive ion etching. 

In this work, test lines of 10-, 30-, 50- (im strip widths 
(W) were fabricated. The length (L) and the spacing to 




Figure 1 Cross-sectional SEM image of typical mesoporous 
silicon. This is fabricated by the anodization of a highly n-doped 
silicon substrate in HF-acetic acid-water solution. 



the surrounding ground plane (S) are set to 1,100 and 15 
(im, respectively. For instance, a structure with 1,100 \im 
in length, 10 \im in width, and 15 \im of spacing to the 
ground plane is denoted as W10S15L1100 as shown in 
Figure 2. 



Results and discussion 

Results 
Electrical test 

These test lines are electrically characterized by the ap- 
plication of a current stress; the experiments have been 
done with a probe station by four wire sensing measure- 
ments. The resistance of a test line is measured versus 
the applied current. The resistance is not constant with 
the current but varies due to temperature dependence. 
Resistance (or resistivity) varies with temperature follow- 
ing a linear law: 



R(T) =£ 0 .(l + TCR 0 .r) 



(1) 



where R(T) is the resistance depending to the 
temperature; R 0 , the resistance at 0°C; and TCR 0 , the 
temperature coefficient of resistance at 0°C. 

Under a current stress, the test line heats, and its re- 
sistance changes, proportional to the temperature rise. 
This resistance variation is linear versus the dissipated 
power (RJ 2 ) as shown in Figure 3. The alpha coefficient 
is defined as the slope and is typical to the Joule heating 
effect. 

For each structure on silicon substrate and porous Si/ 
silicon substrate, the resistance rise is measured, and the 
alpha coefficient is calculated by a linear regression. Ex- 
perimental results are given in Table 1. 




Figure 2 Three-dimensional structure for simulation. Aluminum 
(clear blue) is embedded by the oxide (green) on a silicon substrate 
(dark blue). 



Siegert et al. Nanoscale Research Letters 2012, 7:427 
http://www.nanoscalereslett.eom/content/7/1/427 



Page 3 of 7 



10 



Resistance = 1.1 151 Power + 4.1465 
R 2 = 0.9992 
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Figure 3 The resistance varies with the applied current 
(enclosed). The resistance follows a linear variation with the 
dissipated power. The slope is defined as the alpha coefficient. 



Electrothermal simulation 

The simulator used is the Sentaurus Sdevice tool from 
Synopsis, Inc. (Mountain View, CA, USA). The structure 
is created and meshed. Electrothermal simulation is then 
performed. The model uses Poisson's equation coupled 
with temperature. The temperature induced by the 
current is calculated, and the dissipation is computed by 
using Fourier's law. Both effects are calculated at the 
same time (coupled simulation). 

This model depends on many parameters necessary to 
be implemented in order to make the simulation reliable. 

In order to model the experiments, different parameters 
(dimensional, material properties, boundaries conditions) 
have to be investigated, determined, and implemented in 
the simulation. Electrothermal simulation requires a good 
knowledge of the materials' properties; for this reason, the 
model simulation has to be calibrated. 

• Aluminum parameters. The aluminum parameters 
are defined as follows: the dimensions (width and 
thickness), the resistivity at 0°C, and the TCR at 0°C. 
TCR is determined as 0.0045°C _1 by the following 
measurement. A temperature ramp is applied to 
several patterns (NIST), and their ohmic resistances 
are measured at low current (to avoid joule heating). 
A straight line is obtained, and the TCR is defined as 

Table 1 Alpha coefficient calculation for the different 
dimension strips for both substrates (Si and PoSi) 



Alpha (O W _1 ) 


Dimension strip 




Substrate 






Si 




PoSi/Si 


W10 


1.115 




4.61 


W30 


0.199 




0.93 


W50 


0.09 




0.502 



Table 2 Dimensions (um), resistance (CI), and resistivity 
(CI m) of the three patterns 





Th (um) 


1/1/ (um) 


Mum) 


R (CI) at 30°C 


p (Q m) at 30°C 


W10 


0.94 


9.4 


1,100 


4.053 


3.26 X10" 8 


W30 


0.94 


29.5 


1,100 


1.307 


3.30 X10" 8 


W50 


0.94 


49 


1,100 


0.796 


3.33 X10" 8 



Th, thickness; W, width; L, length. From the resistance (/?) measurement, 
resistivity (p) is calculated. 

the slope divided by the resistance at a given 
temperature. Finally, an average value of TCR is 
obtained. The aluminum thickness is determined by 
means of a focused ion beam (FIB) cut and a SEM 
measurement. The resistivity (p) is calculated by 
taking the resistance measurement at low current 
(with no joule heating influence), the thickness 
measurement, and the theoretical width. Calculation 
is performed for each pattern using Equation 2: 



/?(30°C) = R(30°C) 



T x W 



Results are given in Table 2. The aluminum 
resistivity is almost constant, around 3.30 x 10~* 
m; the resistivity at 0°C is then calculated with 
Equation 3: 



Po 



P(T) 



(1 + TCR 0 ■ T) 



(2) 



(3) 



with T= 30°C by taking the TCR defined as above. 
Stack layer parameters: thermal conductivity and 
thickness. The thermal conductivity coupled with the 
thickness determines the thermal resistance and the 
dissipation. The oxide thickness is determined by the 
FIB cut. Thermal conductivity is taken according to 
the literature and equal to 1.4 W/(K-m)[ll]. 




0.4 
1(A) 

- no calibrated Simulation- • 



0.8 



- Calibrated simulation 



□ Experiments - 

Figure 4 Comparison between experiments (empty squares) 
W10S15L1100 and simulation before (line) and after (broken 
line) calibration. Before the calibration, a mismatch between the 
simulated response and the experiments is found. This mismatch is 
corrected after the determination of the models' parameters. 
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Figure 5 Thermal path for our structure. The main heat path is 
through the substrate. 



Nevertheless, this value can change depending on 
thermal contact resistance between the oxide and 
the substrate. For the silicon substrate, the thermal 
conductivity (X S i) is given by Equation 4 [12]: 



(a + b.T + cT 2 ) 



(4) 



with a = 0.03 K-m-W" 1 , b = 1.56 x 10" 3 m-W" 1 , and 
c=1.65xl0~ 6 m-K^-W" 1 . 

The silicon thermal conductivity (TC) varies with a 
number of parameters such as impurities, doping 
level, and carrier concentration. In this work, we 
have taken the general equation of silicon TC and 
make the assumption that the difference between 
both values was too weak to influence the final 
result (metal heating). 
• Boundary conditions. The meaning of the boundary 
conditions (BC) is to replicate the experimental 
setup [13]. There are two kinds of BCs: thermal and 
electrical. The input of the electrical BC is rather 
trivial; it is just a reference potential, and the input/ 
output terminal (electrical contact) has to be set. 
The specification of the thermal BCs, however, is 
more complex. Two thermal conditions may be 

Table 3 Comparison of the alpha coefficients between 
the experiments and the simulations after calibration 



a Experimental 

(OW _1 ) 



a Simulated 

(OW _1 ) 



Difference (%) 



W10 
W30 
W50 



1.115 
0.199 
0.09 



1.048 
0.180 
0.092 



5.98 
9.33 
2.57 



defined: at the top surface, to take into accounts the 
thermal exchange with air, and at the bottom 
surface, at the chuck/substrate interface. The chuck 
temperature is at 30°C; this temperature is applied 
to the bottom surface. At the top surface, there is 
natural convection with air. We use the traditional 
value of natural convection of 26 W/(K-m 2 ). 

Model simulator calibration on silicon substrate 

An initial simulation is performed with implemented 
values; from the simulation, the resistance versus the 
current are plotted and compared to the experiments as 
shown in Figure 4 (full line). 

Two conclusions can be made. Firstly, values of simu- 
lated resistance and experimental resistances are 
matched at low current. This means that electrical and 
dimensional parameters for the aluminum are correct. 
Secondly, for the high current, the simulation diverges 
from experiments. Simulated DC resistance evolution is 
lower than the experimental DC resistance. The simula- 
tion computes a higher dissipation of the heat generated 
in the line than what actually occurs in the experiments. 

Considering the simplified heat path through the 
structure given in Figure 5, two major parameters can be 
chosen to match the simulation to the experimental 
results. The main path starts from the line through the 
oxide layer down to the silicon and afterward to the 
chuck. Considering the thermal conductivity of silicon 
which is much larger than that of the oxide, the oxide 
thermal conductivity is chosen in order to fit the experi- 
ments. Another parameter, BC at the bottom surface, 
supposed a perfect contact between the chuck and the 
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Figure 6 Thermal path for our structure with a PoSi layer. 
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Figure 7 Electrothermal simulation of W10S15L1 100 on porous silicon/silicon substrate. The porous silicon is simulated by a material layer. 
Thermal conductivity of this layer varies until reaching agreements with experiments. 



silicon substrates. In reality, there are some imperfec- 
tions which existed called contact resistance'. In order 
to take into account this imperfection, a heat coefficient 
is set at this interface. 

Two parameters have been chosen: the thermal con- 
ductivity of the oxide and the BC at the bottom surface. 
The goal is to vary these parameters in the simulation in 
order to reach a fit of experimental data. A design ex- 
periment is performed on these two parameters defined 
as the factors. The response is the alpha coefficient. 
Simulations are performed for each pattern. This makes 
a total of 27 simulations (3x9). 

Finally, having both values unknown, the heat film co- 
efficient and the oxide thermal conductivity are found to 
be 3.4 W/(K-m 2 ) and 1.1 W/(K-m), respectively. We have 
noticed that this last value is not far from the oxide 



thermal conductivity found in the literature, which is 
from 1.0 to 1.2 W/(K m) for a thin silicon dioxide film 
[14]. Simulations were performed for each pattern with 
these values in order to check their validity (Figure 4 
broken line). The comparison of experimental and simu- 
lated alpha coefficients is given in Table 3. 

The parameters found and the conditions in the probe 
station are fixed and used for the PoSi thermal conduct- 
ivity determination. 

PoSi thermal conductivity determination 

With this step, the model simulator is calibrated for our 
experiments (prober station electrical test). The stack 
described previously is now integrated on a 100- (im PoSi 
layer. Therefore, the only unknown is the thermal con- 
ductivity of this new layer in the stack (Figure 6). 




4 6 
Thermal Conductivity SiP (W/(K.m)) 

W10 O W30 O W50 



Exp 



a(W10) = 0.113X 2 -1.750\+9.938 a(W30) = 0.026\ 2 -0.40\+ 2.214 a(W50) = -0.07X+ 0.85 
R 2 = 1.000 R 2 = 1.000 



Figure 8 Alpha coefficient for different thermal conductivities 
(dot) versus the porous silicon thermal conductivity set in 
simulation. From this graph, a porous silicon thermal conductivity 
value is found according to experimental alpha coefficient. 
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Figure 9 Simulations performed for each dimension (W10, 
W30, W50). Simulations fit experiments validating the thermal 
conductivity value determined previously. 
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The thermal conductivity of the porous silicon is largely 
dependent on many parameters, intrinsically to porous 
silicon: pore size, depth, density, substrate P or N [15]. 

Electrical tests are performed on the same test line 
integrated with PoSi. The same methodology is used as 
described above: the resistance versus current is mea- 
sured; the resistance versus dissipated power is plotted; 
and then the alpha coefficient (slope) is computed. 

The value of the mesoporous silicon (our case) ther- 
mal conductivity varies from 1 to 7 W/(K m) in litera- 
ture [16]. In the simulation, a virtual layer is added with 
a variable thermal conductivity (Figure 7); we use three 
levels of variation: 3, 5, and 7 W/(K-m). 

For each thermal conductivity, the alpha coefficient is cal- 
culated. The results for each pattern are given in Figure 8. 

By multi-response optimization, a thermal conductiv- 
ity is deduced to fit the experimental alpha value for 
each pattern. The best results (with the lowest error with 
experimental values) give 4.2 W/(K m). 

To check the determined thermal conductivity value, 
last simulations are performed for each pattern; results 
are given in Figure 9. A good prediction of the resistance 
rise versus the current is found. 



Thermal conductivity micro-Raman measurement 

Micro-Raman scattering spectrometry [17] has been 
used to measure directly thermal conductivity value of 
our porous Si layers. From this measurement, we have 
obtained a value of 5 W/(K m), confirming the validity 
of our simulation methodology. 

The difference between these methodologies can ex- 
plain the small gap between both values. 

Discussion 

In this work, we have evaluated by electrical measure- 
ment coupled with electrothermal simulation the meso- 
porous silicon thermal conductivity. In order to validate 
the determined value, we also performed a micro-Raman 
thermal conductivity measurement. We found a good 
agreement between both values. 

In the model, some assumptions have been taken: 
mesoporous silicon is assumed to be a bulk material. 
The pore size and the morphology are not considered. 
Therefore, the thermal-determined conductivity value is 
an average value. 

Moreover, to apply this methodology, the thermal con- 
ductivity value of the PS needs to be much lower than 
that for the silicon substrate in order to characterize the 
temperature rise difference of the test line between the 
silicon and the porous silicon. The thickness has also to 
be high enough in order to influence the thermal dissi- 
pation and hence the thermal characterization. 



Conclusions 

In this study, we have characterized test lines under a 
probe station by measuring the resistance rise versus the 
applied current. This characterization has been per- 
formed on two different substrates: a silicon substrate 
and a mesoporous silicon/silicon substrate. Firstly, the 
model simulator has been calibrated in order to replicate 
the experiments, and secondly, the thermal conductivity 
of the porous silicon layer has been evaluated. The 
determined value of 4.2 W/(K m) is in agreement with 
the value obtained from micro-Raman scattering mea- 
surements. Our method appears to be efficient and use- 
ful in order to determine thermal conductivity of porous 
silicon morphology provided that its thermal impact is 
significant. 
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